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Abstract. The time evolution of the number density of galaxy clusters and their mass and temperature functions 
are used to constrain cosmological parameters in the spatially flat dark matter models containing hot particles 
(massive neutrino) as well as cold and baryonic matter. We test the modified MDM (A = 0) models with cosmic 
gravitational waves and show that they neither pass the cluster evolution test nor reproduce the observed height 
of the first acoustic peak in AT/T spectrum, and therefore should be ruled out. The models with a non-zero 
cosmological constant are in better agreement with observations. We estimate the free cosmological parameters 
in AMDM with a negligible abundance of gravitational waves, and find that within the parameter ranges h £ 
(0.6,0.7), n £ (0.9,1.1), f v = Q u /Q m £ (0,0.2), (i) the value of Qa is strongly affected by a small fraction of 
hot dark matter: 0.45 < < 0.7 (la CL), and (ii) the redshift evolution of galaxy clusters alone reveals the 
<^^> following explicit relation between J7a and /„: 
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n A + 0.5/„ =0.65 ±0.1. 

This degeneracy is also expected in LSS tests (with a smaller error). The present accuracy of observational 
data allows to bound the fraction of hot matter, /„ £ (0,0.2); the number of massive neutrino species remains 
undelimited, N v = 1, 2,3. 
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1. Introduction by available astrophysical observations. Therefore, assum- 

«2j ■ ing that the initial conditions are given one can derive the 

One of the most important problems of modern cosmology theoretical predictions for LSS formation in a dark matter 

is the formation of large-scale structure in the Universe (DM) model and thcn uge statistical ana i ysis and obs erva- 

(LSS). The last decade progress in observations of LSS tkmal data tQ tegt the viability of the considered models, 

and cosmic microwave background amsotropy (CMBA) . e tQ constrain the allowed range of principa i COS mologi- 

has allowed a productive comparison of theory with obser- ca j p arame ^ ers 
vations. Any realistic cosmological model should be con- 

sistent with LSS observational data in the range of scales Nowadays, the most popular cosmological models are 

from sub-Mpc (Gnedin 1998, Peackock 1997) up to the cold DM wlth n ° n " z ero cosmological constant (ACDM) 

cosmological horizon (Smoot et al. 1992, Jaffe et al. 2001). ( Peebles 1984 ' Kofman & Starobmsky 1985) and mixed 

. ,. t , V i i , , . DM without (MDM) (Fang et al. 1984, Shafi & Stecker 

According to the standard theory of gravitational m- ^ nni Tr , , . . /v. „. ' T . , 

, , it too* i-i tt ■ i i- 1984, Valdarmm & Uonometto 1984, Lukash 1991, 

stability the observable LSS m the Universe has been T . , , . , / a ,,t^-.«-\ 

r , , ,i r ii j • i Lucchm et al. 1996) and with (AMDM (Kahniashvili et 

formed by growth of small density inhomogcneities gen- -,„„„ , T „ . , -,„„„, , V , . \ 

, , . , , , r i-i al. 1996, JNovosyadlyi et al. 1999b and reis. therein) cos- 

erated during the very early stages of cosmological expan- , . , , ^ TT ^n, T ■ • / 

m .„ . , . . i- i . i mological constant, where MDM is m the form oi non- 

sion. Iwo mam features of the primordial perturbation , ° . . ' . . , . 

, , r ,. , ., , baryonic cold (ncutrahno or hypothetical axions) and hot 

held - the Gaussian character of linear density perturba- . . \ . . , 

, , , , toot n j (massive neutrino) particles, 

tions and tne bottom-up LSS formation - are confirmed ' 

The cosmological impact of non-baryonic collisionlcss 

Send offprint requests to: Arhipova N.A. matter depends on the free-streaming path of DM parti- 
e-mail: arna@lukash.asc.rssi.ru cles: at small scales (k > hps) the perturbations smooth 
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out, whereas at large scales (k < kps) they grow to form 
the gravitationally bounded objects of mass greater than 
M F s (Novikov & Zel'dovich 1975). Due to the various free- 
streaming pathes of cold and hot particles (kp S ^> kp S ) 
the growth rate of density perturbations is different in 
CDM and HDM components. The linear perturbation 
power spectrum formed by redshift z, P{k,z), depends 
strongly on the abundance of each component (fi c , f2„) 
and is given by the production of the initial power spec- 
trum, Po(k), the transfer function, T[k, z), and the pertur- 
bation growth rate, D(z) (e.g. Padmanabhan 1993). Both 
the transfer function and the growth factor are DM de- 
pendent and describe the evolution of initial density per- 
turbations during expansion of the Universe (Zakharov 
1979, Bardeen et al. 1986, Holtzman 1989, Eisenstein & 
Hu 1998, Novosyadly et al. 1999a). 

All cosmological models have been re-addressed after 
CMBA experimental detection (Bennett et al. 1996) to 
reveal their positive and negative features. 

Both ACDM and MDM models have met several diffi- 
culties. 

As far as ACDM models are concerned (Kofman et 
al. 1993, Eke et al. 1996, Liddle et al. 1996a,b, Primack & 
Klypin 1996) they demand a high value of the cosmological 
constant > 0.7). In this case ACDM is able to fit 
a set of LSS observational constraints whereas at small 
scales it overproduces the number of collapsed objects by 
a factor of 2 in comparison with the corresponding number 
of gravitationally bounded objects in galactic cataloges 
(Klypin et al. 1996, Gawiser 2000, Bond et al. 2000). 

Regarding MDM models the difficulties are related to 
late galaxy (quasar) formation (Pogosyan & Starobinsky 
1995, Romberg et al. 1996) and too high numbers of ob- 
tained galaxy clusters (Ma 1996, Valdarnini et al. 1998, 
Gardini et al. 1999, Rahman & Shandarin 2001). Standard 
MDM with one, two or three species of massive neutrino 
are ruled out at 2a CL (Novosyadlyj et al. 1999b). 

A way to overcome the sMDM difficulties could be 
the consideration of decaying neutrinos (Bonometto & 
Pierpaoli 1998), or models with cosmic gravitational waves 
(CGW) (Arkhipova et al. 1999, Melchiorri et al. 1999), or 
a non-zero cosmological constant (AMDM) (Kahniashvili 
et al. 1996, Tegmark 1999). In this paper we consider the 
models with stable neutrinos. 

An importance of fundamental CGW has been em- 
phasized by Starobinsky (1979), Rubakov et al. (1982), 
Lukash & Mikheeva (2000). Theoretical predictions for 
cluster abundance in MDM models with CGW has been 
presented by Ma (1996) for the red spectra of density per- 
turbations (n < 1), and by Mikheeva et al.(2001) for both 
blue (n > 1) and red (n < 1) scalar perturbation spectra. 
Although ACDM models with H m < 0.3 normalized by 
COBE 4-year data (Bunn & White 1997) are consistent 
with the cluster number density test, in order to archive 
an agreement with CMBA and cluster abundance data 
in MDM models it is necessary to take into account the 
CGW contribution in the derived value of AT/T at 10° 



angular scale. The latter is estimated by parameter T/S, 
the ratio of the tensor to the scalar mode contributions. 

As an alternative to ACDM and MDM, the AMDM 
models have been considered by Valdarnini et al.(1998), 
Novosyadlyj et al. (1999b), Primack & Gross (2000), 
Andres et al.(2000). The advantage of these models is in 
retaining the inflationary paradigm and the flat Harrison- 
Zel'dovich spectrum (n = 1) with a smaller value (com- 
paring to ACDM) of cosmological A-term. Other notable 
features of AMDM are the possibility of negligible CGW 
contribution (T/S =0) and a small fraction of hot par- 
ticles: even 10% of massive neutrino in matter content 
(fl v /il m ~ 0.1) could change the value of the cosmologi- 
cal constant, being however in good agreement with other 
independent tests: CMB data (Hu et al. 2000), distant 
SNIa (Perlmutter et al. 1999), Hipparcos data (Feast & 
Catchpole 1997), QSO lensing (Kochanek 1996). 

In this paper we consider MDM (with CGW) and 
AMDM (with T/S=0) applying to these models the clus- 
ter abundance and evolution tests. The evolution test con- 
strains most efficiently the parameter f2 m , and hence 
in spatially flat cosmological models (e.g. Bahcall et al. 
1997, Donahue &Voit 1999, Henry 2000). Our aim is to 
demonstrate how the presence of hot matter influences 
the estimation of Da, h, and other cosmological parame- 
ters. Here, instead of doing an exact \ 2 analysis, we rather 
look for the tendencies and correlations between cosmo- 
logical parameters when introducing a fraction of massive 
neutrino in the Universe. 

We describe our models in Section 2, and the galaxy 
cluster tests in Section 3. The results are discussed in 
Section 4, with the conclusions in Section 5. 

2. Cosmological models with massive neutrino 

We assume that DM is given by mixture of CDM and 
HDM components in the flat background space. The free 
model parameters are: 

— f2 m , the total matter density in the Universe 
(fl m = 1 — ft\ = fift + O c + Cl v , the latter being 
density parameters of baryons, cold, and hot particles, 
respectively); 

— f v = Q^/Qjtj, the fraction of hot DM; 

— N v , the number of massive neutrino species; 

— h , the Hubble constant in units 100 km s _1 Mpc _1 ; 

— n , the slope-index of post-inflationary density pertur- 
bation power spectrum. 

The massive and massless neutrinos^] are described by 
the corresponding distribution functions which are eval- 
uated from the Boltzmann-Vlasov collisionless equations. 

1 All neutrinos are active, the total number of neutrino 
species is three. Accordingly, N v = 1, 2, 3. 
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Cold particles are treated as a pressureless fluid (pc = 0) . 
Baryons and photons are described as an ideal hydrody- 
namic fluid satisfying the Euler equations of motion. Here 
we choose the fixed value for baryon density parameter, 
f2& = 0.015//i 2 . All DM components interact with each 
other only gravitationally. 

The calculation of the perturbation dynamics (the 
transfer functions) demands a joint solution of the de- 
scribed self-consistent set of equations. It can be done nu- 
merically (Seljak & Zaldarriaga 1996) or analytically in 
long and short wave regions (Zakharov 1979, Bardeen et 
al. 1986, Holtzman 1989), the possibility of semi-analitical 
approximations of some transfer functions has been shown 
by Eisenstein & Hu (1998), Novosyadlyj et al.(1999a). 

Assuming the power low post-inflationary density per- 
turbation spectrum, Po(fc)oofc™, the final power spectrum 
of total density perturbations can be expressed as 



P(k,z) = Ak n T\k,z) 



D{z) 



D(0) 



(1) 



where A is the normalization constant, T(k, z) is the total 
transfer function, D(z) is the growth factoip], 



D(z) = 



g(n m (z)) 
l + z 



(2) 



and g(Q m (z)) is the suppression coefficient (Kofman & 
Starobinsky 1985). 

According to Carroll et al.(1992) the function g(x) can 
be approximated as 



g(x) — 5x \2x 7 + 



2 + .x(209 -x)\ 1 



70 



(3) 



abd the current matter abundance Q m (z) is written as 
follows: 

<W*) - ^ x-nj+Vr^ ■ n ™ s a " (0) ' (4) 

We use the transfer function approximations for sCDM 
(Bardeen et al. 1986) and AMDM (Eisenstein & Hu 1998) 
models. 



3. Cluster mass and temperature functions and 
the evolution test 

The mass function for the gravitationally bounded halos 
of mass greater than M formed in the flat Universe by 
redshift z is given by (Press & Schechter 1974) 



N(> M, z) = / n(M , z)dM 



(5) 



M 



2 The function D(z) is given as the growth rate of linear 
density perturbations in the model without hot matter (with 
the same parameters fl m , h, and n). 



where n(M, z)dM is the comoving number density of col- 
lapsed objects with masses lying in the interval (M, M + 
dM): 

n(M,z) = Jl*L J l*^V=4g . (6) 
1 ' ' V 7T M a 2 (R,z) 1 dM 1 V ' 

M = f Trpi? 3 , p is the mean matter density, and S c is the 
critical density contrast for a linear overdensity able to 
collaps. The rms amplitude of density fluctuation in the 
spheres of radius R at redshift z is related to the power 
spectrum as 



a 2 {R, z) 



1 

2^2 



P(k,z)\W(kR)\ 2 k 2 dk 



(7) 



where W(kR) is a Fourier component of the top-hat win- 
dow function: W(x) = -^(sinx — xcosx). We denote 
a R = o{R, 0) for z = 0. 

In the matter-dominated Universe S c = 1.686 (Eke et 
al. 1996, Viana & Liddle 1996), in the flat models with A- 
term S c depends weakly on the current matter abundance 
(Liddle et al. 1996a, Locas & Hoffman 2000). 

The theoretical mass functions obtained with the help 
of the Press-Schechter formalism (eqs.(5,6,7)) arc in good 
agreement with other methods including numerical simu- 
lations (Efstathiou et al. 1988, Lacey & Cole 1994, Gross 
et al. 1998, Eke et al. 1998, Bode et al. 2000, Wu 2000). 
Due to the exponential dependence of n(M, 0) on an (see 
eq.(6)) the cluster number density is very sensitive to the 
value of as • 

Rich galaxy clusters are strong X-rays sources charac- 
terised by the gas temperature. The physical mechanisms 
of clusters formation promp the relation between cluster 
mass and temperature, T s ooM 2 / 3 , confirmed by numeri- 
cal simulations (e.g. Navarro et al. 1995). The proportion- 
ality coefficient depends slightly on cosmological model 
(mainly, on the parameter O m ). Here, we use the T — Al 
relation for isothermal gas given by Eke et al. (1996): 



7.75 



6. 



13 V5A" + 3 



Cl m (z) 



^] (l + z)Kev(8) 



where /3(~ 1) is the ratio of the galaxy kinetic energy to 
the gas thermal energy, AT(~ 0.76) is the hydrogen mass 
fraction, and the mass M15 is given in units 1O 15 /i _1 M0. 
The value A cr is the ratio of a mean halo density (within 
the virial radius of collapsed object) to the critical density 
of the Universe at corresponding redshift. For Q, m < 1, A cr 
can be derived analytically and approximated as A cr = 
178f2^ 45 for z = 0. 

As well as the cluster mass function N(> M) = N(> 
M, 0), the X-cluster temperature function N(> T) = N(> 
T, 0) is also sensitive to an. The comparison between 
the observed temperature (Henry & Arnaud 1991) and 
mass (Bahcall & Cen 1993) functions shows a good agree- 
ment with eq.(8) at z = 0. Comparing with observational 
data both tests provide a powerful constraint on a- value 
in different DM models (Valdarnini et al. 1998, Bahcall 
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& Fan 1998, Durrer & Novosyadlyj 2001). On the other 
hand CMBA strongly depends on the density perturba- 
tion spectrum as well. The consistency with all these tests 
determines the model parameters (e.g. Bridle et al. 1999, 
Mikheeva et al. 2001). 

The cluster mass (temperature) function and its evo- 
lution in ACDM has been discussed in detail (e.g. Viana 
& Liddle 1998, Bode et al. 2000, Henry 2000). To achieve 
a better consistency with cluster evolution data the value 
of n m should be lower than that obtained from the cluster 
number density test at z = 0. 

The dependence of the cluster mass (temperature) 
function on z originates due to different growth rate of 
density perturbations in DM models. It is reflected in the 
appearance of the suppression factor g(Q m (z)) in models 
with Q\ ^ (see eqs.(3,4)): 



a(R,z) 



ctr g{Vt m (z)) 
1 + z g(Q m ) 



(9) 



Below, we consider the cosmological models containing 
massive neutrino. In this case z-dependence is also present 
in the transfer function due to the free-streaming of hot 
matter particles (the effective growth rate becomes scale 
dependent). For the mass estimation of rich galaxy clus- 
ters (see eq.(6)), the effect of free-streaming is negligible. 

4. Results and Discussion 

4.1. MDM models with non-zero tensor mode 

To test MDM models with a zero cosmological constant 
we follow the normalization procedure given by Arkhipova 
et al. (1999), Mikheeva et al. (2001). All models are nor- 
malized by as by the best fit present day cluster mass 
(Bahcall & Cen 1993, Bahcall et al. 1997) and tempera- 
ture (Henry & Arnaud 1991, Donahue & Voit 1999 ) func- 
tions. To achieve an agreement with AT/T data (Bennet 
et al. 1996) the contribution of cosmic gravitational waves 
is required. The value of the derived parameter T/S de- 
pends mainly on the spectral index n, the abundance of 
hot matter Sl v , and the Hubble constant h. 

Fig.l shows the present-day cluster temperature func- 
tions N{> T) for different values of Cl v in MDM models 
with CGW normalized by as ~ 0.52. The normalization 
does not depend on f2„, therefore all curves cross each 
other at some fixed point corresponding to the mass M in 
the sphere of radius 8/i _1 Mpc. 

A significant contribution of both CGW and massive 
neutrinos is needed to fit the large scale CMBA observa- 
tional data. The consideration of the red spectra (n < 1) 
cannot reduce T/S substantially. E.g., in MDM models 
with fixed Q v = 0.2, h = 0.6, and n = 0.9,1,1.1 : 
T/S=0.6,1.7,3.6 respectively, which strongly suppresses 
the height of the first acoustic peak in AT/T spectrum 
(Mikheeva et al. 2001). 

In this respect, higher values of f2„ would be preferable 
for MDM models: an increase of f2„ > 0.2 reduces the pa- 
rameter T/S, however the problem arises with small-scale 



I i i i | i i i | i i i— 
MDM 

<7 a normalization: a^O^Q^ - 52 * 0130 ^ 




0.4 



0.6 0.8 
T [KeV] 



Fig. 1. The present day cluster abundance N(> T) 
in MDM models normalized by as = 0.52 with Q, v = 
0,0.1,0.2,0.3 (solid, dot, dash, long-dash lines, resp.), 
Q, b = 0.015//i 2 , h = 0.6, n = l. The needed T/S is shown. 
The data points by Henry & Arnaud (1991). 



clustering and the Ly a -cloud formation tests (Gnedin 
1998). On the contrary, at low f2„ < 0.2 neither possible 
changes of h nor models with three species of massive neu- 
trino can decrease sufficiently a high contribution of CGW 
into CMBA at large angular scales. The latter obviously 
leads to a low height of the first acoustic peak inconsistent 
with the BOOMERANG observations (Netterfield et al. 
2001). Fig. 2 demonstrates that the introduction of even 
a small cosmological constant would drastically improve 
the situation with MDM-CGW model parameters. 

But let us return again to MDM models without A- 
term. Fig. 3 presents the evolution of galaxy clusters. None 
of the considered models can fit the data at high z > 0.3. 
We varied h (Fig. 3a), (Fig. 3b), and n (not shown). 
The evolution tracks are practically independent of the 
post-inflationary spectral index for n € (0.9, 1.1). Even ne- 
glecting the observational problems with Ly^-forest, the 
case with high Q, w > 0.2 cannot help to achieve an agree- 
ment with the evolution of galaxy clusters in MDM domi- 
nated models. All the models predict the number of galaxy 
clusters at high z > 0.3 at least two orders of magnitude 
smaller than the observed one. This fact indicates that the 
evolution of density perturbations should be slower than 
that found in MDM without a cosmological constant. 
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Fig. 2. The introduction of a cosmological con- 
stant in MDM-CGW models normalized by erg = 
0.52^7°- 52+013Om) with n v = 0.1, fl b = 0.015/fr 2 , h = 
0.6: (a) T/S as a function of n for fl A = 0,0.3,0.4,0.5 
(solid, dot, dash, long-dash lines, resp.), and (b) the height 
of the first acoustic peak, i?i = 70 ^ , as a function of 
r^A for n — 0.9,0.95,1,1.05,1.1 (solid, dot, dashed, long- 
dashed, dashed-dot lines, resp.). 

4.2. MDM models with non-zero cosmological 
constant 

Let us consider flat AMDM models without CGW nor- 
malized by the COBE 4-year data. This allows an ex- 
plicit derivation of crs(cmb) as a function of model param- 



Fig. 3. The cluster evolution N(> M = 8 • 10 14 M Q ,z) 
in MDM models normalized by cr% — 0.52 with n = 1, (a) 
f2„ = 0.2, h = 0.5,0.6,0.7 (solid, dot, dash lines, resp.), 
and (b) h = 0.6, Q v = 0,0.1,0.2,0.3 (solid, dot, dash, 
long-dash lines, resp.). The needed T/S is shown. The data 
points by Bahcall & Fan (1998). 



eters. On the other hand, from observational abundance of 
galaxy clusters it is possible to obtain (cl) as a function 
of n m , e.g. cr 8 (cl) = 0.52r2-°- 52+013n "« (Eke et al. 1996, 
Liddle et al. 1996b). The coincidence between as(cmb) and 
as (cl) would limit the parameter 17 a by ~ 0.7 in ACDM 
(Bahcall & Fan 1998, Henry 2000) and - 0.5 - 0.6 in 
AMDM (Valdarnini et al 1998, Novosyadlyj et al. 1999b) 
models. 
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Fig. 4. The cluster evolution N(> M = 8- 10 14 M o 
ACDM models normalized by ag = 0.52fir. 
with n A = 0.7, n b = 0.015/ft 2 , (a) n = 1, h = 
0.65,0.70,0.75 (dot, dash, long-dash lines, resp.), and (b) 
h = 0.7, n = 0.8,0.9,1.0,1.1 (solid, dot, dash, long-dash 
lines, resp.). The data points by Bahcall & Fan (1998). 



According to Henry (2000) the cluster evolution test 
has a stronger upper-limited parameter f2 m in ACDM 
models than does the cluster abundance test at z — 0. As 
an example we demonstrate the evolution of the number 
density of galaxy clusters of mass M > 8 • 10 14 M Q within 
the comoving radius 1.5/i _1 Mpc in ACDM model with 
£1a = 0.7, as functions of h (Fig. 4a) and n (Fig. 4b). The 
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COBE normalization 



(b) 

n H /n M =o.2 
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n A =o.3i 
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h=0.6 n=1.0 



_L 



0.4 



0.6 0.8 
T [KeV] 



Fig. 5. The present day cluster abundance A^(> T) in 
AMDM models normalized by COBE 4-year data with 
f^A = 0,0.31,0.45,0.65 ( dot-dash, long-dash, dash, dot 
lines, resp.), il;, = 0.015//i 2 , h — 0.6, n = 1, and /„ = 0.1 
(a), 0.2 (b). The data points by Henry & Arnaud (1991). 

cluster evolution is not practically affected by changes of 
the spectral index. Regarding the parameter h, the models 
with h > 0.6 are preferable. 

While the requested value of cosmological constant in 
ACDM models is quite high, a small amount of hot parti- 
cles in AMDM models leads to the agreement with obser- 
vations for a smaller Q\. 

We performed the calculation of functions A^(> T) for 
different parameters J1a and f v (Fig. 5). As we can see, 
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to achieve an agreement with today cluster temperature 
function, we need a lower value of cosmological constant 
in the models with a higher fraction of hot matter. For 
/„ G (0,0.2) ,h G (0.6,0.7) ,n G (0.9,1.1) , the value of 
cosmological constant remains in the range ( la CL): 



0.35 <ft A < 0.7 



(10) 



Taking into account the Ly a forest test, models with f v ~ 
0.1 are preferable. In these case parameter tljy is limited by 
0.45 < fl A < 0.65 (in agreement with some other results, 
cf. Valdarnini et al.(1998), Novosyadlyj et al.( 1999b), 
Primack & Gross (2000), Durrer & Novosyadlyj (2001)). 

To obtain the redshift evolution functions of galaxy 
cluster abundance, the normalization of the power spec- 
trum is done according to the present-day cluster con- 
centration, as (cl), which depends on Q\. At la level the 
best fit for as (cl) is not sensitive to the hot particle abun- 
dance, h G (0.6,0.7), and n G (0.9,1.1). Figs. 6,7 present 
the galaxy cluster evolution in AMDM models as function 
of parameters J7a, fui and h. Changing h and n (within 
their ranges) does not practically influence the evolution 
tracks. Increasing f v reduces the needed value of £l\ (in 
agreement with the cluster number density test). 

For h G (0.6,0.7), n G (0.9,1.1), the value of n A de- 
pends mainly on the fraction of hot matter /„. We find 
the following approximation between the parameters £l\ 
and /„ G (0, 0.2) from the cluster redshift evolution alone 
(la CL): 



A + 0.5.^ = 0.65 ±0.1 



(11) 



In connection with this degeneracy two points should be 
emphasized. 

* As we see, the redshift evolution test slightly ex- 
aggerates the value of in comparioson with the 
present day cluster abundance (cf. eq.(10)). To be 
in the agreement with both tests the value of the 
cosmological constant in AMDM models should be 
found as n A G (0.45,0.7). 

* Eq.(ll) can also be understood as a proportionality 
between the values of VL m and f v : more matter leads 
to a greater fraction of hot matter and vice verce. Such 
a degeneracy is physically clear and could be straight- 
forwardly observed in LSS tests. E.g., Novosyadlyj' et 
al. (1999b) indicate the correlation between Vl m and 
Vt v from the cluster power spectrum and Ly a data. 
We guess that a joint LSS analysis will confirm eq.(ll) 
with a better precision. 

In conclusion, we stress that the present observational 
data constrain only the parameter f u G (0, 0.2), the num- 
ber of species of massive neutrino (JV„ = 1,2,3 ) remains 
undetermined. 

5. Conclusions 

Our main conclusions are as follows: 
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Fig. 6. The cluster evolution A^(> M = 8 • 10 14 M Q , z) in 
AMDM models normalized by a$ = 0.52rjL°- 52+013am) 
with fl A = 0,0.31,0.45,0.55,0.65 (solid, dot, dash, dot- 
dash, long-dash lines, resp.), Sit, = 0.015//i 2 , h — 0.6, n = 
1, and ft, — 0.1 (a), 0.2 (b). The data points by Bahcall 
& Fan (1998). 



The introduction of massive neutrino in spatially flat 
dark matter models affects drastically the estimation 
of the cosmological parameters from galaxy cluster 
observational data. 

MDM models with cosmic gravitational waves and 
a negligible A-term do not pass the observational 
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The introduction of even a small cosmological con- 
stant in MDM-CGW models essentially improves the 
situation with the model parameters, indicating that 
A-term is a more powerful instrument (than T/S) to 
reconcile theory with observations. 

The cluster abundance and/or evolution tests alone 
clearly hint at the existence of a non-zero cosmological 
constant in a set of spatially flat AMDM models with 
a negligible amount of CGW. A rise of T/S cannot 
cancel the A-term. However, the opposite is not true: 
the estimation of T/S in AMDM cosmologies requires 
better data and a more subtle analysis. 

The presence of a small fraction of massive neutrino, 
f u £ (0, 0.2) , (i) reduces the required value of cosmo- 
logical constant (in comparison with ACDM), (ii) is 
preferable from the point of view of the ly a forest test, 
and (iii) does not discriminate the estimate of other 
cosmological parameters, h £ (0.6,0.7), n £ (0.9, 1, 1), 
N v = 1,2,3. 
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Fig. 7. The cluster evolution N(> M = 8 • 10 14 M Q ,z) 
in AMDM models normalized by er 8 = 0.52ft m a52+013nm 
with n A = 0.55, n b = 0.015/h 2 , n = l,(&)h = 0.6, f v = 
0,0.1,0.2 (solid, dot, dash lines, resp), and (b) /„ = 0.1, 
h = 0.6,0.65,0.7 (dot, dash, long-dash lines, resp.). The 
data points by Bahcall & Fan (1998). 



- For f v £ (0,0.2), h € (0.6,0.7), n £ (0.9,1.1), the 
value of cosmological constant remains in la ranges: 

* f^A £ (0.35,0.7) (from nearby cluster abundance), 

* f^A £ (0.45,0.75) (from cluster redshift evolution), 

* fl A £ (0.45,0.7) (from the both tests). 

— We find the following approximation between param- 
eters J7a and f v £ (0, 0.2) from the cluster redshift 
evolution alone: 



0.5/„ = 0.65 ±0.1 , (lcrCL). 



— If massive neutrinos constitute only ~ 10% of the to- 
tal DM density (f^ = 0.1) the models with fl\ £ 
(0.5, 0.65) satisfy the observational data best. 

It may well be that the DM nature is more com- 
plex than just AMDM with a negligible amount of CGW. 
Moreover, even in the framework of this simple model 
we cannot delimit today the fraction and the number of 
species of massive neutrino, the accuracy of the available 
data is still low to do it. Nevertheless, we can assert that 
the abundance of cosmological massive neutrinos is small, 
f u < 0.2, corresponding less than few eV to the neutrino 
mass. Also, today we can conclude that 

* the value of the estimated cosmological constant is 
very sensitive to even a small fraction of massive 
neutrinos, and 



constraints: they require a significal contribution of 
both (i) CGW (T/S >1) and (ii) hot matter (/„>0.2). 
The first condition reduces the first acoustic peak 
in AT/T spectrum by a level inconsistent with the 
BOOMERANG measurements, the second condition 
contradicts the Ly Q -system formation test. 



* the continuing progress in LSS cosmological observa- 
tions assists in solving the DM problem. 
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